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•^ ', Abstract. We established a new eighth-order iterative method, con- 

f^^ I sisting of three steps, for solving nonlinear equations. Per iteration the 

method requires four evaluations (three function evaluations and one 
evaluation of the first derivative). Convergence analysis shows that this 
-^ [ method is eighth-order convergent which is also substantiated through 

^ I the numerical works. Computational results ascertain that our method 

r^ ' is efficient and demonstrate almost better performance as compared to 

the other well known eighth-order methods. 
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1. Introduction 



This paper concerns the numerical solution of non-linear equations 
of the general form f{x) = 0. Such equations appear in real world 
situations frequently while there is no closed form solution for them. 
That is why the numerical solution of these types of equations draw 
c^ ■ much attention to itself day by day. One of the common problems 

encountered in science and engineering problems is that given a single 
variable function f{x), find the values of x for which f{x) = 0. The 
root of such nonlinear equations may be real or complex. There are 
two general types of methods available to find the roots of algebraic 
and transcendental equations. First, direct methods, which are not al- 
ways applicable to find the roots, and second, iterative methods based 
on the concept of successive approximations. In this case, the general 
procedure is to start with one or more initial approximation (s) to the 
root and attain a sequence of iterates, which in the limit converges 
to the true solution. Multipoint iterative methods for solving nonlin- 
ear equations are of great practical importance since they overcome 
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theoretical limits of one-point methods concerning the convergence or- 
der and computational efficiency. Here, we focus on the simple root 
of nonlinear scalar equations by iterative method. Let the function 
/ : D C 3fJ — > 3fJ be a sufficiently differentiable function and a G -D be 
a simple root of f{x) = 0. The famous Newton's method [4j of order 
two which can be defined as x^+i = x„ — /(a;„)//'(x„) is one of the 
oldest and the most applicable method in the literature. Traub [11] 
proposed the concept of efficiency index as a measure for comparing 
methods. This index is prescribed by p^^"", where p is the order of 
convergence and n is the whole number of evaluations per iteration. 
Kung and Traub [9j then presented a hypothesis on the optimality of 
roots by giving 2""^ as the optimal order. This means that the Newton 
iteration by two evaluations per iterations is optimal with 1.414 as the 
efficiency index. By taking into account the optimality concept many 
authors have tried to build iterative methods of optimal higher order 
of convergence. Anyway, these schemes are divided into two main cat- 
egories. First the derivative free and second, the methods in which 
(ffist, second, ...) derivative evaluation is used per cycle. In this paper 
we present a new eighth-order iterative method to find a simple root a 
of the nonlinear equation f{x) = 0. We will compare our new method 
with well known existing eighth-order methods, namely proposed in 
[20], [12], [15], [H, [22], [8], [7], [3], [21], . 

This paper is organized as follows: In section 2, we describe the 
eighth-order method and prove that the method obtained preserves 
their convergence order. This new method agree with the Kung and 
Traub conjecture for n = 4. In section 3, we will briefly state the well 
known established methods in order to compare the effectiveness of the 
new method. Finally, in section 4, the new method is compared in the 
performance with some well known eighth-order methods. Numerical 
results indicate that the our proposed method give better performance. 
Consequently, we have found that the new eighth-order method is con- 
sistent, stable and convergent. 



2. Development of the method and analysis of 

convergence 

In this section, we will define a new eighth-order method. In order 
to establish the order of convergence of this new method, we state fol- 
lowing definitions: 
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Definition 2.1. Let f(x) be a real function with a simple root a and let 
Xn be a sequence of real numbers that converge towards a. The order 
of convergence m is given by 

lim ;^"+^~" =C^0, (2.1) 

where ( is the asymptotic error constant and m G -R"*". 



Definition 2.2. Let (5 be the number of function evaluations of the new 
method. The efficiency of the new method is measured by the concept 
of efficiency index [121 [TU] and defined as 

^'l'. (2.2) 

where \i is the order of the method. 



2.1. New Eighth-order Method. In this section, we construct our 
eight-order method by considering three step cycles. Kung [16] devel- 
oped a one-parameter family of fourth-order methods, which is written 
as 

Vn — Xn -J-, r, l^-'JJ 

/ (a^n) 
^ _, /(x„) + /3/(y„) / /(y„) \ 

" ^" f{Xn)^{(i-2)f{yn)\r{Xn))' ^'-^^ 

where /3 is a constant. In particular, the special method for /3 = —1/2 
is as follows: 

J \^n) /OCX 

f \Xn) 

_ 2f{xn)-f{yn) f f{yn) \ , . 

" ^" 2f{xn)-hf{yn)\r{xn))- ^'-'^ 

To achieve an S*'' order method, we use the first and second steps of 

the three step cycle from the above method and in the third step we 

apply the Newtons iteration 

Xn+l = Zn 777 T- (2.7) 

As we can see, this method consists of three evaluation of the func- 
tion and two evaluation of the first derivative per iteration, which has 
efficiency index 1.5157. To derive a scheme with a higher efficiency in- 
dex, we approximate f'{zn) using a Hermite interpolation polynomial 
which is equal to the function f{x) in the domain D [D is the interval 
in which / has a simple root). 
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To approximate f'{zn), we construct a Hermite interpolation poly- 
nomial, H{x), that meets the interpolation conditions 

H{Xn) = f{Xn), H{yn) = fiVn), H {Zn) = f {Zn) and H'{Xn) = f'{Xn)- 

The H{x) could be defined as follows: 

H{x) = UJo{x)f{Xn) + UJl{x)f{yn) + UJ2{x)f{Zn) + UJo{x)f'{Xn), (2.8) 

where the interpolation basis functions ujo{x), uji{x), ijJ2{x) and cZ;o(a;) 
are cubic polynomial function that satisfy in the following conditions 

(i) UJo{Xn) = l^UJoiVn) = 0,Wo(^n) = 0,Wo(^n) = 0? 
(ii) Ui{Xn) = 0,Ui{yn) = l,Ui{Zn) = 0,u[{Xn) = 0, 

(m) uj2ix.a) = 0,W2(l/n) = 0,a;2(^n) = l,W2(a;„) = 0, 

(iv) UJo{Xn) = l,OJo{yn) = 0,CJo(^n) = 0,UJQ{Xn) = 1. (2.9) 

From (?), we have uJo{x) = A{x — yn){x — Zn){x — B), where A, B are 
constants. The conditions uJo{xn) = 1 and ojQ^Xn) = imply 

t ^Xn Zn yn jj \Xn yn)\Xn ^nj 

[Xn yn) \Xn ^nj '^X-n Z^, |/n 

Thus 



yx Xn)\X Zfi) 

UJo{x) 



\Xn ynjyXn Zn) 



yx Xnjy^Xn Zfi Vn) 



(2.10) 



\Xn ynjyXn Zn) 

From iii), we have uji{x) = C{x — x„)^(a; — Zn) where C is a constant. 
The condition a;i(|/„) = 1 implies 

1 



Thus 



Similarly, 



C 



U;,{x) = ; ", "\ . (2.11) 



(l/n- 


XufiVn - 


ZnY 


{x 


Xn) [X 


~ Zn) 


iVn 


-XnY{yn 


~ Zn) 


(x 


Xn) \X 


-Vn) 



^2{x) = ;_^'':,y -\ . (2.12) 

\Zn Xn) \Zn Vn) 



From (if), we have uj{x) = D{x — Xn){x — yn){x — Zn) where D a 

g 
1 



constant. The conditions ujQ{xn) = 1 gives 



D 



\Xn yn)\Xn Zn) 

Hence 



— / \ l-^ Xn)\X yn)[X Zn) , . 

ujq{x) = — — ^ . (2.13) 

\Xn yn)\Xn Zn) 
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With the specific expressions of the interpolation basis functions u)o{x) 
uji{x) uj2{x), aJo(x), we have 



Hf,)- i^-yn){x-Zn) 



\Xn VnjyXn Zn) 



yX XfijyZXji yn Zj^ 
yXfi yn)\Xn Zn) 



fM 



\X Xn) \X Zj^J yX Xn) yX yn) r / 

yVn XnlyUn ^nj \Zn X^j yZj^ y^j 

\X XnjyX ynjyX Znj ^ , , 
H 7 w ~ ^ J [Xn)i 

\Xn yn)\Xn ^nj 



(2.14) 



Differentiating the above equation then putting x = Zn and Simplifying 
we can get 



Xn Zn yn Z^ X^ yn 

- Zn) f fiVn) - fiXn) \ _ (x„ - Zn) 
{yn - Xn) \ {yn " Xn) j {Vn " Xn) ' 



, \Xn Zn) I J yVn) J \Xn) I \Xn Zn) „i , n „// n 

+ 7 T \ —( ^ r ~ 7 \i (^") ~ / (^«)' 



(2.15) 
and subsequently, we can find 

-n \Zn) -^I [Xni Zn\~r J [ynj Zn\ J [Xn, yn\~T~ \Xn Zn) J [yn, Xn, Xn\ / \X ) , 

(2.16) 
where f[xn, Zn] (similarly f[yn, Zn] and /[x„, yn]) and /[y„, x„, x„] are 
defined below by taking into consideration the divided differences: 

f{Xn) - f{Zn) 



J [Xn, Zn 
J [yn, Xn, Xn 



Xn Zn 
f[yn,Xn] - f'{Xn) 



yn Xn 

Now replacing f'{zn) in the equation (12.71) by H'{zn) from the equation 
f l236|) . we find that 

Xn+l ^ J [Xn, Zn\ ~\~ J [l/n; ^nj J [Xn, yn\ T \Xn Zn) J [yn, Xn, Xn\ J \X), 

(2.17) 
Now this method consists of three evaluation of the functions and one 
evaluation of the first derivative per iteration, now it has improved 
efficiency index 1.6817. Now we prove the convergence of this method 
by the following theorem: 

Theorem 2.1. Let us consider a as the simple root of the nonlinear 
equation f(x)=0 in the domain D and assume that f(x) is sufficiently 
smooth in the neighborhood of the root. Then, the iterative scheme 
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defined by i \2.5L (12 .Gp and ( I2.17P is of local order eight and has the 
following error equation 

Cn+i = {44 - clc3C4,)el + 0{el), 
where e„ = a;„ — a and Ch = ^j , h=l,2,3.... 

Proof. We provide the Taylor series expansion of each term involved in 
(12.51) . (12. 6 p and (I2.17p . By Taylor expansion around the simple root in 
the n^^ iteration, we have 

f{.Xn) = /'(a)[e„ + C2el + c^el + c^e^ + clel 

+ ceel + cjel + cge^ + Cc^el + cioe^° (2.18) 

+ ciie^^ + cise^^ + 0(e^3)] 
and, we have 

f'{xn) = f{a)[l + 2c2e„ + Scge^ + Ac^el + hc\el 

+ Gcge^ + 7c7e^ + Scge^ + Ocge^ + lOcioe^ (2.19) 

+ llcnei° + 12ci2e;,i + 0(e;,2)]. 
Further more it can be easily find 

ii^ = e„ - c^el + 1c<^el - 2c,el + + 0(e^2). (2.20) 

By considering this relation and equation (12. Sp . we obtain 

yn = a + C2el + 2(c3 - cl)el + + 0(6^^). (2.21) 

At this time, we should expand /(?/„) around the root by taking into 
consideration (12.210 . Accordingly, we have 

fiVn) = f'{a)[c2el + 2(-c^ + cs)el + + 0(e^^)]. (2.22) 

Using ([MID, (EIISD, fl^^ and (1^:^ in the equation ([2SD, we can 
obtain 

z^ = a- C2Cset + {^4 + 2clc^ - Ic^c^ - 2cl)el + .... + Oie";^). (2.23) 

On the other hand, we have 

/(z„) = f[a)[-C2C^ei + (^-clc^-2c,C2-2cl)el + + 0(e^^)]. (2.24) 

Now we expand the Taylor series of each existing divided differences in 
the denominator of equation (I2.17p . We obtain 

/K, Zn] = f'{a)[l + C2en + csel + c^el + + 0{e^^)l (2.25) 

f[yn,zn] = f{a)[l + C2el + (2C2C3 " 2c^)e^ + + 0{e^^)], (2.26) 



(2.28) 
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f[xn,yn] = f'ia)[l+C2en+icl+C3)el+{-2cl+3c3C2+Ci)el+ +0{el% 

(2.27) 

f[yn, Xn, Xn] = /'(«) [Cg + 2C3e„ + (3C4 + €203)6^ 

+ (-2c^C3 + 2C4C2 + 24 + Ac5)el + .... + 0{el^)]. 

By consider the above mentioned relations (I2.25p -( l2.28p in the equation 
(I2.17p . we can find 

3 

Xn+l = (C2C3 - C2C3C4)e^ + [-3C2C3 + -C2C4 - 4C2C3 + 8C2C3C4 

+ iclcl - 2c5clc3 - 2clcl - Ac2clc4\e^ + ... + 0{e\^). 

D 

This imphes that the order of convergence for this method is eight i.e. 
optimal order of convergence and its efficiency index is 8^^"^ ~ 1.6817, 
which is more that 1.4142 of Newton's method, and 1.5650 of three- 
step methods |S1 12], and is equal to 1.6817 of [20], [12], US], HI, [22], 
|H], [Z], [3], |2I] but numerical performance is better to almost all the 
eighth-order methods ( later shown in the Tables 2-8). 

3. Well established eighth-order Methods 

First we are giving here some well established eighth-order methods: 

3.1. Bi et al. Methods, [ [20], 2009]: 
Method I 

y. -"" 



Xri. 



f'iXn) 



Vn- 



Method II 

^n i/n ~ 
X-n+l ^n' 

Method II 

^n 1/n 



2/(x„) - f{yn) 

V{Xn) - 5/(l/„) 

/(x„) + (7 + 2)/(z„) 



f{yn) 



fM +'lf{Zn 



f{Zn 



J Y^ri) yn\ ~r J [Zny Xn, Xn\\Zn yn) 

(3.1 



2/(y„) f f{yn) 



2/(a;„) \fixn 
/(x„) + (7 + 2)/(z„) 

f{Xn)+7f{Zn) 



V{yn) ffiyn) 



+ 



fiVn) 
J [Xn/ 



f{yn) 



f{Zn 



J i^nj yn\ ~r / [^^n, -^n, Xn\ \Zn yr, 

(3.2) 



2/(x„) \f{Xn) 



f{yn) 

f{Xn) 



f{yn) 

nxn) 
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•^n+1 ^n 



/(x„) + (7 + 2)/(z„) 



f{Xn)+7fiZn) 



f{Zn) 



Method IV 



(3.3) 



Zn Vn 



2f{xn) - 3f{y„] 



•^n+l ^n 



f{Xn 

f(Xn) + {7 + 2)fiZny 
f{.Xn)+lf{Zn) 



-2/3 



fiVn) 
f{Zn 



(3.4) 
where 7 G i? and denominator is not equal to zero. 



3.2. Sharma et al. Methods, [ [12], 2010]: 

Method I 



Zn Vri 



f{Xr. 



Method II 



[f{xn)-2f{y„) 
2 



1 + -ir, — r + 7 



fiXn) 



Zn = y 



f{Xn] 
f{Xn) 



Xn+1 Zn 

Method III 



L/(xO-2/(i/„) 

/(x.) + (7 + l)/(Zn) 
fM + lf{Zn) 

f{Xn) 



fiVn 
I'M. 

f[Xn,yn]fiZn) 
J [yni Zn\ J [Xm Zn\ 

fiVn 



Z-n Vn 



f'iXr. 
f[Xn,yn]fiZn) 
J [yni Zn\J [Xn-i Zn 

fiVn) 



[f{Xn)-2f{yn)\ \f'{Xn) 



Xn+1 Zn 



1 + 7 



fjZn 
f{Xr. 



^"'^ f f[Xn.yn]f{Zn) 
J [yn 1 Zn\J [Xn , Zn 



where 7 G i? and denominator is not equal to zero. 



(3.5) 



(3.6) 



(3.7) 



3.3. Thukral Method, [ [15], 2010]: 

f{Xn? + f{yn? 



^n •^n 



f'iXn)ifiXn)-fiyn)y 



Xn+1 2^?! 



1 + yu: 



2\ 2 



IJ'Tl 



2{^^n?~Q{^^nf+^^"''^ ' ^^^'"^ 



+ 4- 



f{yn) f{Xn) 



f{Zn) 
f'iXn] 
(3.8) 
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where //„ = {f{yn)/f{xn)) and denominator is not equal to zero. 



3.4. Wang et al., [ [2l], 2010]: 

^n Vn 



fiVn) 



•^n+l ^n 



2f[Xn,yn\ - f'{Xr. 



(3.9) 



3.5. Sragolzaei et al. Method, [ [H], 2011]: 



^n i/n 



-d + i^ 



■'^n+l ^n 



f{Zn) 



'2f[Zn, Xn] + f[Zn, yn] " 2/[|/„, Xn] + (yn - Zn)f[yn, Xn, Xn] ' 

(3.10) 



3.6. Cordero et al. Method, [ [3J, 2011]. 



Un Zn 



Zn Xfi 



f{Zn) 
I'M 



fjXn 

nxr. 



f{Xn) - f{yn) 

f{xn)-2f{yn) 



f{Xn) - fiVn) , 1 f{Zn) 



/(x„)-2/(y„) 2/(y„)-2/(z„)J 

Un Zn 



Xn+l — Un <J „^ 



f'iXn) 



Vn Xn 



(3.11) 



3.7. Soleymani Methods, [([7], [5]); (2011, 2012)]: 

Method I 

\_^f'{yn?-f'{Xn)'] fiXn) 



Zn Xn 



Xn-\-l Zn 



Method II 



f'iyny J f'iXn) 

fjZn) 

J (.yn) ~r -^J [Znt Xn, Xn\\Zn yn, 

2f{Xn) 



(3.12) 



Zn Xn 



"^n Zn 



Xn-\-l "^n fu 



f'{Xn) + f'{yny 
f{Zn) 



f'iyn) 



J'{Xn){:if'(yn)-f'(Xn)) 

f'{xn)+f'{y„) 



(3.13) 
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3.8. Kim Method, [ [22], 2012]: 



Vn 



l + {(3-2)Un + fiulJ f'iXn) 



^ ^ _ /' l + aUn + bVn \ fjZn) .^ ^^. 

\l + CUn + dVnJ f'{Xn) + f[yn,Xn,Zn]{Zn-Xn)' 

where «„ = /(l/„)//(x„) , Vn = f{zn)/ f{xn) , /3, A, /i, a, 6, c, and d 
are constant parameters and related by /3 = (A — /x — 2/3) /2, a = —2, 
c = —3, (i = 6 — 3. 

4. Numerical Testing 

To demonstrate the performance of the new eighth-order method, 
we take seven particular non-linear equations. We wiU determine the 
consistency and stability of results by examining the convergence of 
the new iterative method. We will give estimates of the approximate 
solution produced be the eighth-order method. 

Here we consider, the following test functions to illustrate the accu- 
racy of new iterative method. The root of each nonlinear test function is 
also listed in from of each page up to fifteen decimal places, when such 
roots are non-integers. All the computations reported here we have 
done using Mathematica 8, Scientific computations in many branches 
of science and technology demand very high precision degree of numer- 
ical precision. The test non-linear functions are listed in Table- 1. The 

Table 1. Test functions and their roots. 

Non-linear function Roots 

/i(a;) =sin(x) - ^fo 

Ux) = ^-x'-j, + l 1 

f^{x) = exp(sm(x)) - 1 - f 

/4(x) = X + sin(^ ) 

f,(x) = v/^^T8sin(;^) + -fi^ - v^+ A .2 

U{x) = cos(x) - X 0.739085133215160 

fr{x) = exp(x) + cos(x) -1.7461395304080124 

results of comparison for the test function are provided in the Table 2- 
8. It can be seen that the resulting method from our class are accurate 
and efficient in terms of number of accurate decimal places to find the 
roots after some iterations. 
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Table 2. Errors Occurring in the estimates of the root 
of function /i by the method described with initial guess 
xo = 0.7. 



Methods 




fiixi) 


fi{x2) 


fiixs) 


(2.17) 




0.695e-5 


0.654e-60 


0.336e-665 


(13. ip with 7 = 


1 


0.791e-3 


0.363e-18 


0.740e-95 


(I3.2P with 7 = 


1 


0.205e-3 


0.428e-21 


0.169e-109 


(33D with 7 = 


1 


0.628e-3 


0.115e-18 


0.233e-97 


(MP with 7 = 


1 


0.887e-2 


0.177e-8 


0.140e-28 


(3.5P with ^ = 


1 


0.689e-4 


0.339e-48 


0.139e-535 


(3.6P with -/ = 


1 


0.821e-4 


0.233e-47 


0.224e-526 


(3.7P with ^ = 


1 


0.754e-4 


0.918e-48 


0.790e-531 


(31D 




0.129e-2 


0.174e-28 


0.257e-261 


(XTU) 




0.275e-4 


0.239e-53 


0.511e-593 


(3.14) with (A 


/i,6) = 


= (0,0,4) O.lOle-3 


0.414e-46 


0.201e-512 


(XTH) 




0.148e-4 


0.116e-65 


0.512e-860 


(3.12) 




0.133e-2 


0.601e-10 


0.548e-32 


(3.11) 




0.830e-4 


0.409e-39 


0.702e-357 


dsD 




O.lOle-4 


0.414e-58 


0.217e-645 
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Table 3. Errors Occurring in the estimates of the root 
of function /2 by the method described with initial guess 
xo = 1.2. 



Methods 



f2{Xl) f2{x2) f2{x3) 



(2.17|) 




O.llle-3 


0.497e-31 


0.823e-250 


(KTh with 7 = 


1 


0.372e-3 


0.709e-24 


0.638e-161 


(1S2D with 7 = 


1 


dgt 


dgt 


dgt 


(13.3) with 7 = 


1 


0.165e+0 


.136e-3 


0.442e-15 


(13.4) with 7 = 


1 


0.768e-2 


0.116e-6 


0.427e-21 


(13.5) with 7 = 


1 


0.317e-2 


0.192e-18 


0.376e-148 


(13.6) with 7 = 


1 


0.547e-2 


0.141e-16 


0.308e-133 


fl3.7) with 7 = 


1 


0.420e-2 


0.176e-17 


0.186e-140 


dM) 




0.213e+l 


0.536e-l 


0.271e-7 


(3.10]) 




0.109e-l 


0.431e-14 


0.362e-113 


f3.14|) with (A 


fi,b) = 


= (0,0,4) 0.496e-4 


O.lOle-33 


0.318e-271 


(3.13]) 




0.667e-3 


0.119e-24 


0.128e-198 


(13.12^ 




0.545e-4 


0.419e-9 


0.247e-19 


(13.11^ 




0.992e-2 


0.683e-15 


0.467e-120 


(ED 




0.713e-3 


0.409e-24 


0.494e-194 
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Table 4. Errors Occurring in the estimates of the root 
of function f^ by the method described with initial guess 
Xq = —0.55. 



Methods 



/3(a;i) f3{x2) MX3) 



(EUD 

( IXTD with 7 = 1 
(1S2D with 7=1 
f l33|) with 7 = 1 
(^ with 7 = 1 
f l33|) with 7 = 1 
f l33|) with 7 = 1 
( K7\f with 7 = 1 



(ETOD 

f lXTij) with (A,/i,6) 

(ESD 
(EH 
dSTI]) 



0.628e-2 
0.175e-l 
0.266e+l 
0.176e+0 
0.486e-2 
0.520e-2 
0.883e-2 
0.684e-2 
0.271e+l 
0.107e-l 
(0,0,4) 0.623e-2 
0.737e-2 
0.424e-2 
0.176e-l 
O.llOe-2 



0.344e-20 

0.682e-9 

0.260e+2 

0.736e-4 

0.243e-5 

.212e-18 

0.141e-16 

0.186e-17 

0.626e+3 

0.344e-15 

0.337e-18 

0.568e-17 

0.198e-5 

0.284e-14 

0.141e-24 



0.168e-184 

0.549e-146 

0.527e+3 

0.871e-17 

0.614e-12 

0.173e-149 

0.652e-135 

0.610e-142 

0.225e+5 

0.582e-123 

0.261e-148 

0.774e-138 

0.431e-12 

0.162e-116 

0.104e-199 
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Table 5. Errors Occurring in the estimates of the root 
of function f^ by the method described with initial guess 
xo = 0.1. 



Methods 




h{xi) 


h{x2) 


U{x,) 


(2.17|) 




0.467e-14 


0.371e-147 


0.370e-1478 


dSH) with 7 = 




0.340e-12 


0.287e-115 


0.629e-1043 


(K7\f with 7 = 




0.199e-10 


0.174e-97 


0.513e-881 


(13. 3 p with 7 = 




0.124e-ll 


0.128e-109 


0.173e-991 


(13. 4 p with 7 = 




0.274e-4 


0.568e-15 


0.508e-47 


(13. 5 p with 7 = 




0.756e-ll 


0.106e-91 


0.156e-738 


(13. 6 p with 7 = 




0.756e-ll 


0.106e-91 


0.162e-738 


(IMP with 7 = 




0.756e-ll 


0.106e-91 


0.159e-738 


(ESD 




0.141e-9 


0.322e-80 


0.242e-645 


(3.10P 




0.354e-10 


0.130e-85 


0.435e-689 


(3.14P with (A 


fi,b) = 


= (0,0,4) 0.169e-ll 


0.145e-97 


0.442e-786 


(3.13P 




O.lOle-10 


0.146e-90 


0.283e-729 


(3I2D 




0.472e-3 


0.999e-8 


0.447e-17 


(anD 




0.681e-ll 


0.423e-92 


0.950e-742 


(ESD 




0.256e-ll 


0.618e-96 


0.706e-773 
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Table 6. Errors Occurring in the estimates of the root 
of function /s by the method described with initial guess 
Xq = -3. 



Methods 




Mxi) 


h{x2) 


hix-s) 


(2.17) 




0.283e-3 


0.555e-26 


0.108e-207 


(13. ip with 7 = 


1 


0.756e-2 


0.218e-7 


0.129e-29 


(I3.2P with 7 = 


1 


0.690e-2 


0.143e-7 


0.243e-30 


(33D with 7 = 


1 


0.745e-2 


0.201e-7 


0.107e-29 


(3.4p with 7 = 


1 


0.179e-l 


0.159e-3 


O.lOle-7 


(3.5P with ^ = 


1 


0.333e-3 


0.599e-25 


0.659e-199 


(3.6P with -/ = 


1 


0.321e-3 


0.447e-25 


0.633e-200 


(3.7P with ^ = 


1 


0.328e-3 


0.519e-25 


0.209e-199 


(31D 




0.251e-3 


0.409e-25 


0.199e-199 


(XTU) 




0.380e-3 


0.157e-24 


0.138e-195 


(3.14) with (A 


/i,6) = 


= (0,0,4) 0.432e-3 


0.545e-25 


0.344e-200 


(XTH) 




0.705e-3 


0.362e-23 


0.173e-185 


(3.12) 




0.622e-2 


0.676e-5 


0.779e-ll 


(3.11) 




0.122e-2 


0.507e-23 


0.229e-185 


dsD 




0.307e-3 


0.132e-25 


0.160e-204 
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Table 7. Errors Occurring in the estimates of the root 
of function /g by the method described with initial guess 
Xq = 1.5. 



Methods 




/6(Xl) 


f(i{x2) 


hix-s) 


(2.17) 




0.696e-6 


0.176e-55 


0.300e-452 


(13. ip with 7 = 


1 


0.682e-2 


0.124e-10 


0.132e-45 


(I3.2P with 7 = 


1 


0.706e-2 


0.141e-10 


0.227e-45 


(33D with 7 = 


1 


0.686e-2 


0.126e-10 


0.143e-45 


(MP with 7 = 


1 


0.650e-l 


0.383e-3 


0.129e-7 


(3.5P with ^ = 


1 


0.415e-5 


0.990e-48 


0.103e-388 


(3.6P with -/ = 


1 


0.421e-5 


O.llle-47 


0.261e-388 


(3.7P with ^ = 


1 


0.418e-5 


0.109e-47 


0.164e-388 


(31D 




0.249e-4 


0.203e-40 


0.385e-329 


(XTU) 




0.142e-5 


0.222e-51 


0.797e-418 


(3.14) with (A 


/i,6) = 


= (0,0,4) 0.228e-5 


0.317e-50 


0.454e-409 


(XTH) 




0.234e-5 


0.185e-50 


0.289e-411 


(3.12) 




0.129e-l 


0.311e-5 


0.180e-12 


(3.11) 




0.893e-5 


0.142e-45 


0.589e-372 


dsD 




0.870e-6 


0.363e-54 


0.332e-441 
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Table 8. Errors Occurring in the estimates of the root 
of function fj by the method described with initial guess 
Xq = —2.3. 



Methods 




/7(Xl) 


f7{x2) 


/7(X3) 


(2.17) 




0.563e-6 


0.167e-54 


O.lOle-442 


(13. ip with 7 = 


1 


0.105e-2 


0.464e-14 


0.179e-59 


(I3.2P with 7 = 


1 


0.529e-2 


0.301e-ll 


0.318e-48 


(33D with 7 = 


1 


0.490e-3 


0.222e-15 


0.946e-65 


(3.4p with 7 = 


1 


0.578e-2 


0.554e-6 


0.507e-14 


(3.5P with ^ = 


1 


0.412e-4 


0.749e-39 


0.893e-317 


(3.6P with -/ = 


1 


0.443e-4 


0.135e-38 


0.977e-315 


(3.7P with ^ = 


1 


0.428e-4 


O.lOle-38 


0.960e-316 


(31D 




0.718e-3 


0.478e-28 


0.186e-229 


(XTU) 




0.788e-4 


0.774e-37 


0.671e-301 


(3.14) with (A 


/i,6) = 


= (0,0,4) 0.231e-4 


0.551e-42 


0.580e-343 


(XTH) 




0.788e-6 


0.357e-53 


0.647e-432 


(3.12) 




0.384e-3 


0.269e-8 


0.132e-18 


(3.11) 




0.246e-4 


0.910e-41 


0.314e-332 


dsD 




0.592e-5 


0.357e-46 


0.630e-376 
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5. Conclusion 

In this work, we have developed a new eighth-order convergent method 
for solving non-linear equations. Convergence analysis shows that our 
new method is eighth-order convergent which is also supported by the 
numerical works. This iterative method require evaluation of three 
functions and one first derivative during each iterative step. Compu- 
tational results demonstrate that the iterative method is efficient and 
exhibit better performance as compared with other well known eighth- 
order methods. 
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